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INTRODUCTION 

Recently  there  have  been  numerous  studies  [1)-[61  on  the  dielectric 

loaded  gyrotron  for  a  wide  band  application.  It  has  been  found  [2] ,  [4] ,  [6] 

that  there  exist  three  unstable  inodes  characterized  by  their  axial  phase 

velocities  (vph);  the  long  wavelength  mode  (LWM,  vp^  >  c) ,  the  intermediate 

-1/2 

wavelength  mode  (IWM,  c  >  vp^  >  ce  ),  and  the  short  wavelength  mode 

(SWM,  v^  <  c).  Here  e  and  c  are  the  dielectric  constant  and  the  velocity  of 
light.  For  a  small  axial  velocity  spread  (<  1%),  two  slow  wave  modes  (IWM 
and  SWM)  yield  very  promising  bandwidth  capability  [2],  especially  when  two 
modes  are  mixed  by  placing  the  beam  close  to  the  axis  [4] ,  16] .  However,  the 
nature  of  the  SWM  [2] ,  [4] ,  whose  perturbed  fields  are  almost  entirely 
supported  by  the  electron  beam,  not  by  the  waveguide,  raises  the  difficulties 
related  to  the  excitation  and  collection  of  the  electromagnetic  waves 
[2],  [4].  One  possible  solution  to  this  difficulty  is  to  utilize  the  dielectric 
material  as  a  center  rod  [7],  rather  than  as  an  outside  wall  loading  [  1] - [5] . 

By  using  dielectric  material  as  a  center  rod,  it  will  be  shown  that  the 
frequency  range  of  the  IWM  is  extended,  while  the  contribution  of  the 
troublesome  SWM  is  minimized. 

The  dispersion  relation  for  the  gyrotron  with  a  dielectric  center  rod 
has  been  derived  by  authors  [7]  for  general  azimuthal  (6)  harmonic  numbers 
(i.e.  I  $  0).  In  this  paper,  we  will  perform  a  detailed  numerical  investigation 
of  this  dispersion  equation  for  a  broad  range  of  the  physical  parameters,  and 
compare  the  results  with  those  of  the  gyrotron  with  the  outer  dielectric 
loading  12],  (4],  [6].  For  simplicity,  the  present  investigation  is  limited 
to  the  azimuthally  symmetric  (i.e.  I  -  0) ,  transverse  electric  (TE)  perturbations. 
Moreover,  in  view  of  the  relative  insensitivity  of  the  fast  wave  mode  (LWM) 
to  the  system  parameters  12],  16],  17],  we  will  devote  our  attention  only  to 
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the  slow  wave  nodes  (IWM  and  SWM) .  A  complete  parametric  optimization 
process  for  wide  bandwidth  will  be  carried  out,  assuming  that  the  axial 
velocity  spread  of  the  beam  electrons  is  small  (=  1%).  The  optimization 
is  carried  out  in  the  spirit  of  maximizing  the  contribution  of  the  IWM  to 
the  bandwidth,  and  minimizing  that  of  the  SWM.  In  addition,  the  perturbed 
field  profiles  are  examined  in  order  to  distinguish  the  IWM  and  SWM. 

The  slow  wave  modes  (IWM  and  SWM)  for  the  center  rod  gyrotron  do  not 
compete  each  other  in  their  contribution  to  the  instability  as  much  as  for 
the  wall  clad  configuration.  This  cooperative  nature  of  the  IWM  and  the 
SWM  results  in  an  extended  IWM  region  compared  to  the  wall  clad  geometry. 

That  is,  the  SWM  begins  to  contribute  significantly  at  higher  frequency  than 
it  does  in  the  outside  loaded  gyrotron.  On  the  other  hand,  it  will  be  shown 
that  the  bandwidth  itself  for  the  center  rod  configuration  is  approximately 
as  wide  as  that  for  the  wall  clad  one  (see  Sec.  III).  Thus  by  putting  the 
dielectric  rod  at  the  center,  we  are  able  to  maintain  about  the  same  wide  band 
capability,  while  eliminating  some  of  the  difficulties  associated  with  the  SWM. 

A  brief  review  for  the  derivation  procedure  of  the  dispersion  relation 
will  be  given  in  Sec.  II.  The  expressions  for  the  perturbed  fields  are  also 
given  in  Sec.  II  for  later  use.  The  optimization  for  the  wide  bandwidth  with 
1%  of  the  axial  velocity  spread  is  carried  out  in  Sec.  III.  The  physical 
parameters  to  be  optimized  are  the  thickness  ratio  of  the  dielectric  center  rod 
(Rd/Rc)»  the  dielectric  constant  (e),  the  conducting  wall  radius  (Rc) ,  and  the 
beam  center  location  (RQ) .  In  Sec.  IV,  the  perturbed  field  profiles  are 
investigated  in  order  to  examine  the  individual  contribution  of  the  IWM  and  the 
SWM  to  the  bandwidth.  Especially  we  compare  these  field  profiles  with  those 
of  the  wall  clad  configuration.  The  summary  of  the  comparison  is  given  in  the 
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DISPERSION  RELATION 

The  cross  section  of  the  gy rotron  with  a  dielectric  center  rod  is 
shown  in  Fig.  1.  A  cylindrical  dielectric  (e)  rod  of  radius  R^  is  located 
concentrically  with  the  conducting  wall  of  radius  Rc«  The  hollow  electron 
beam  passes  through  the  space  between  the  dielectric  center  rod  and  the 
conducting  wall.  The  individual  electrons  undergo  the  cyclotron  motion 
with  Lamor  radii  r^  about  the  beam  center  location  Rq,  under  the  influence 
of  the  constant  applied  axial  magnetic  field  Bq.  In  addition,  the  electron 
beam  moves  downstream  with  the  axial  velocity.  The  cylindrical  coordinates 
(r,  6,  z)  are  employed. 

The  dispersion  relation  is  derived  within  the  framework  of  the 
Vlasov-Maxwell  equations  for  the  fields  EOc.t)  and  B(x,t),  and  for  the  beam 
electron  distribution  function  f(x,£,t).  Here  x,  £,  and  t  refer  to  the 
spatial,  momentum,  and  the  time  coordinates.  Further,  any  quantity  ip  is 
linearized  according  to 

«(x,t)  *  *Q(r)  +  ^(r)  expfifkz  -  «t)J  (1) 

with  the  equilibrium  quantity  ip^  and  the  small  Fourier  decomposed  perturbation 
♦j.  Note  that  we  limit  our  attention  to  the  azimuthally  symmetric  perturbation 
(3/39=0)  with  the  frequency  w  and  the  axial  wavenumber  k.  Moreover,  we  will 
consider  the  transverse  electric  (TE,  Ezl  =  0)  perturbation  only.  It  is 
assumed  that  the  beam  is  tenuous,  and  the  beam  thickness  is  small  [1],  [2], 

14] ,  16] ,  [7] .  In  order  to  examine  the  effect  of  the  axial  velocity  spread, 
the  equilibrium  distribution  function  is  assumed  to  be  Lorentzian  [2] ,  [4] , 
[6] ,  [7]  in  the  axial  momentum  pz>  that  is, 

f0  *  P2  AKPZ  _  Pz>2  +  Pz  A2!’1.  (2) 
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Here  pz  is  the  average  axial  linear  momentum  and  A  is  the  axial  momentum 

2 

spread  ratio.  The  beam  is  further  assumed  to  be  monoenergetic  with  ymc  ,  and 
the  average  transverse  (axial)  velocity  is  given  by  (cBz) • 

Since  the  details  of  the  procedures  in  obtaining  the  dispersion  relation 
are  given  in  Ref.  7,  here  we  present  only  the  outline.  Making  use  of  the 
thin  beam  approximation  and  the  boundary  conditions  on  the  azimuthal  electric 
field  Eqi  at  r*0,  R^,  R^,  and  R^,  within  a  normalizing  factor  we  obtain  the 
perturbed  fields  and  . 


'61 


JjCy), 


2  lAJN  Vx>  ~  AJJ  Ni^x^* 

B 

-  \  INj(xc)  J1(x)-J1(x£.)  NjCx)]. 
3 


0  $  r  <  R 


R.<r  <  R„ 
d  -  0 


Rn<  r<R 
O'  ~  c 


(3) 


*rl  -  -  £  E.l-  Bil 


r  7  !?  <*.!>■ 


(4) 


The  jump  condition  on  Bzl  across  the  beam  furnishes  the  desired  dispersion 
relation.  Namely, 

2  2 

B  ve;  c 

JL  a  -  -  ’ 

BD  2y  Rq2  ♦  ijkj  cBz  yA/y^]2 

with  the  wave  admittance  11] ,  12] ,  [4] ,  [6] ,  (7]  BN/Bp  is  given  by 
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bn  =  2B1'  bd  = 


-  **0  B2  B3’ 


B1  =  W  AJN  '  N1(XC)  AJJ, 


B2  =  W  AJN  "  AJJ* 


(6) 


B3  =  W  N1^X<P  '  W  VV>> 


AJJ  -  ,  „ 

Jl(xd) 

w 

A  E  yd  J0(yd^ 

ajn 

NjCXj) 

■Xd 

W 

W 


The  arguments  of  the  first  (J)  and  the  second  (N)  kind  Bessel  functions  are 


(7) 


In  Eq.  (5),  the  Doppler-shifted  beam  mode  u>B  is  defined  by 

u>„  =  kcB.  ♦  us  h, 

D  1  C 


(8) 
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2  2 

the  Budker  parameter  (v)  is  given  by  v  h  Ne  /me  ,  u>c  =  eBQ/mc  is  the  non- 

2  -1/2 

relativistic  electron  cyclotron  frequency,  and  y  =  ( 1  -  6  )  is  the 

z  z 

axial  mass  factor.  Here  N  is  the  total  number  of  electrons  per  unit  axial 
length,  and  (-e)  and  m  are  the  charge  and  the  rest  mass  of  the  electron. 

It  is  easy  to  show  that  the  dispersion  relation  in  Ref.  7  yields  the  identical 
result  (Eq.  (5))  in  the  limit  of  t= 0. 

The  dispersion  relation  (5)  is  numerically  solved  for  the  axial  wave- 
number  k  in  terms  of  the  frequency  w,  when  other  beam  and  geometric  parameters 
(v,  ,  6Z>  A,  Rq,  c,  Rc»  Rj)  are  given.  If  the  solution  k  is  complex,  the 

mode  is  unstable,  and  the  gain  is  given  by  the  negative  value  of  the 
imarginary  part  of  k  (i.e.  -k^  =  -  Im(k) ,  when  k^  <  0).  The  characteristics 
of  the  unstable  modes  thus  obtained  are  very  similar  to  those  with  the  wall- 
clad  dielectric  gyrotron  [2].  The  TE  perturbation  exhibits  three  unstable 
modes;  one  fast  wave  mode  (LWM)  and  two  coexisting  slow  wave  modes  (1WM  and 
SWM)  separated  by  a  stable  band  near  u=ck  line.  Both  the  LWM  and  IWM  are 
originated  from  the  unstable  coupling  of  the  beam  mode  w  (Eq.  (8))  and 
the  beam- free  waveguide  mode  wG  (the  solution  of  =  0  in  Eqs.  (5)-(6)). 

On  the  other  hand,  the  SWM  is  driven  by  a  quasi-static  magnetic  dipole 
moment  instability  [6],  [8],  characterized  by  highly  localized  perturbed 
fields  near  the  beam  location  ([2],  [4],  [6]  and  see  Sec.  IV).  As  in  the 
case  of  the  wall  clad  geometry,  the  reduction  in  the  gain  due  to  the  axial 
velocity  spread  (A  in  Eq.  (2))  is  the  least  for  the  LWM,  moderate  for  the  IWM, 
and  the  largest  for  the  SWM.  On  the  other  hand,  the  bandwidth  at  small  axial 
velocity  spread  (A  <  1%)  is  the  broadest  for  the  SWM,  intermediate  for  the 
IWM,  and  the  narrowest  for  the  LWM.  Since  our  objective  is  to  achieve  the 
wide  bandwidth,  we  will  concentrate  only  the  slow  waves  (IWM  and  SWM)  in  the 
remainder  of  this  paper.  In  view  of  the  difficulties  associated  with  the 
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SWM  (excitation  and  collection  of  the  electromagnetic  waves) ,  we  attempt  to 
find  the  parameter  conditions  where  the  IWM  shows  a  wide  bandwidth,  over 
which  the  SWM  is  substantially  suppressed. 

In  light  of  our  intention  to  minimize  the  SWM  contribution,  there  are 
several  ways  to  distinguish  the  integrated  two  slow  waves  (IWM  and  SWM). 

One  method  is  to  utilize  the  different  vulnerability  of  their  gains  on  the 
velocity  spread  (A) .  Although  we  choose  A=l%  for  our  investigation,  we 
therefore  examine  the  gain  for  A=3%  as  well.  If  the  gain  is  substantial  for 
both  A=l%  and  3%,  we  identify  this  instability  due  to  the  IWM.  On  the  other 
hand,  if  the  gain  is  greatly  reduced  for  A=3%,  we  label  them  as  the  SWM. 

The  other  method  to  distinguish  the  two  modes  is  to  examine  the  perturbed 
field  profile.  If  the  field  profile  is  very  similar  to  that  of  the  beam-free 
waveguide,  then  the  instability  is  due  to  the  beam-waveguide  coupling  IWM. 

If,  however,  the  field  profile  is  highly  localized  near  the  beam  location, 
we  attribute  the  gain  to  the  SWM.  The  former  method  is  used  in  Sec.  Ill  and 
the  latter  in  Sec.  IV. 


13/14 


NSWC  TR  81-387 


PARAMETRIC  OPTIMIZATION 

In  this  section  we  will  obtain  the  optimized  physical  parameters  for 
the  wide  band  gyrotron  amplifier;  the  thickness  ratio  of  the  dielectric 
center  rod  (Rj/R,.) >  the  conducting  wall  radius  (Rc) ,  the  dielectric  constant 
(e),  and  the  beam  location  (Rq).  We  again  emphasize  that  the  optimization 
is  for  the  bandwidth  due  to  the  IWM,  minimizing  the  contribution  of  the 
troublesome  SWM.  In  the  remainder  of  this  paper,  we  assume  the  following 
beam  parameters. 

6,  =  0.4,  0  =  0.2,  v  =  0.002,  (9) 

corresponding  to  the  60.3  KV  of  the  anode  voltage  and  6.8  Amp.  of  the  total 
axial  current.  For  future  reference,  we  also  define 

R°  =  4.197  c/uc.  (10) 

It  can  be  shown  that  R®  is  the  optimized  wall  radius  when  the  dielectric  is 
absent  [2] ,  [4] ,  (6] .  The  axial  velocity  spread  (A)  of  the  beam  is  assumed 
to  be  small  (A=l%).  However,  in  order  to  examine  the  contribution  of  the  SWM 
(to  be  minimized),  case  for  A=3%  is  also  examined  as  for  reference. 

Since  the  IWM  results  from  the  coupling  of  the  beam  mode  Wg(Eq.  (8)) 
and  the  beam- free  waveguide  mode  (BN=0  *n  Ec*'  much  informati°n  on  the 

IWM  can  be  extracted  by  examining  the  two  modes,  n)g(61  ,Bz,u>c)  and  »RC»0 • 

In  Fig.  2,  a  schematic  diagram  of  the  mode  characteristics  is  given  in  the 
space  of  the  conducting  wall  radius  (R^)  and  the  dielectric  constant  (e).  The 
dispersion  curves  for  the  parameters  denoted  by  A,  B  and  C  are  schematically 
shown  in  the  upper  comer.  When  the  beam  mode  Ug  (the  straight  line  in  the 
upper  comer  diagram)  grazes  the  waveguide  mode  (case  B) ,  the  IWM  gain  is 
maximum.  For  the  parameters  above  the  maximum  gain  curve  (solid  curve  in  Fig.  2), 
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nig  intersects  nore  than  one  time  (case  C) ,  and  for  the  parameters  below 
the  maximum  gain  line,  uig  does  not  intersect  at  all  (case  A).  Since  the 
INM  gain  is  maximum  at  the  u>g-u>G  intersecting  points,  it  is  obvious  that  the 
SUM  is  dominant  when  the  parameters  (R^  and  e)  are  well  below  the  maximum 
INM  gain  line,  while  at  near  and  above  it,  the  INM  is  dominant  as  shown. 

Above  the  maximum  INM  gain  line,  where  tog  intersects  oiG  more  than  one  time, 
the  gain  yields  multiple  maxima  in  w-space  with  a  valley  in  between.  Therefore, 
we  expect  that  the  maximum  bandwidth  line  (broken)  is  located  slightly  above 
the  maximum  gain  line  (solid).  Along  the  maximum  bandwidth  line,  the  gain 
at  the  valley  in  k^-u>  diagram  is  just  high  enough  so  that  the  bandwidth  covers  the 
both  maxima,  to  yield  broadest  bandwidth.  As  one  moves  further  above  the 
maximum  bandwidth  line,  the  valley  of  k^-u  diagram  is  too  deep  to  extend  the 
bandwidth  to  both  maxima,  thereby  abruptly  decreasing  the  bandwidth.  All 
these  arguments  will  be  later  confirmed  in  the  numerical  investigations. 

After  seeing  that  the  beam-waveguide  grazing  condition  (case  B  in  Fig.  2) 
plays  an  important  role  in  predicting  both  the  gain  and  the  bandwidth  of  the 
INM,  we  now  proceed  to  find  the  grazing  conditions  for  the  beam  parameters 
given  in  Eq.  (9).  The  results  are  summarized  in  Fig.  3,  where  the  values  of 
the  wall  radius  (R£,  solid  lines)  and  the  dielectric  constant  (c,  broken 
lines)  that  make  the  beam  mode  graze  the  waveguide  mode  (upper  comer) 
are  plotted.  The  thickness  ratio  of  the  dielectric  rod  (R^RC)  varies  from 
0.10  to  0.35  with  0.05  increment.  The  results  in  Fig.  3  are  similar  to  those 
for  the  wall  clad  configuration  (Fig.  7  in  Ref.  2) .  However,  there  are  several 
noteworthy  differences.  In  general,  the  center  rod  configuration  yields 
larger  wall  radius,  and  smaller  dielectric  constant  compared  to  the  wall  clad 
gyrotron.  For  example,  for  the  same  thickness  ratio  (R,j/Rc  *  0.20  for  the 
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center  rod,  Rw/R£.  =  0.80  for  the  wall  clad,  where  is  the  inner  radius  of  the 
dielectric  material),  R./R*?  =  0.7  -  0.8,  e  *  12  -  13  in  the  IWM  frequency 

C  C 

range  for  the  center  rod,  and  R /R°  *  0.5  -  0.6,  e  *  13  -  15  for  the  wall 

c  c 

clad  case.  The  larger  wall  radius  for  the  center  rod  configuration  nay  be  an 

advantage  when  the  operating  frequency  is  very  high.  We  note  from  Fig.  3 

that  for  the  thickness  ratio  Rd/R£  *  0.20,  the  grazing  Rc  or  e  remains 

relatively  unchanged  over  the  IWM  frequency  range  (i.e.  u/uc  >  1.1).  That  is, 

at  this  dielectric  configuration  (Rd/Rc  *  0.20),  the  beam  mode  grazes  or 

nearly  grazes  the  waveguide  mode  for  a  broad  range  of  the  frequency  in  the 

IWM  region,  thereby  resulting  in  a  wide  bandwidth.  We  therefore  conclude 

that  the  optimized  rod  thickness  ratio  is  given  by  Rd/Rc  *  0.20.  However, 

we  also  note  that  the  value  of  e  or  Rc>  although  relatively  flat,  varies 

more  sensitively  than  that  for  the  wall  clad  case  [2).  This  enhanced 

sensitivity  aids  the  center  rod  gyrotron  in  achieving  wider  bandwidth, 

although  at  smaller  gain,  compared  to  the  wall  clad  case.  This  will  be 

confirmed  later  in  the  numerical  investigations. 

The  relative  sensitivity  of  the  gain  to  the  variation  of  R£  or  e  for 

the  center  rod  configuration  necessitates  additional  optimization 

investigations  on  R£  and  e.  These  are  illustrated  in  Figs.  4  and  5.  With 

the  beam  parameters  (9)  and  the  optimized  rod  thickness  ratio  R./R-  ■  0.20, 

0  c 

the  linear  gain  (-k^)  is  numerically  obtained  from  the  dispersion  relation 
(5),  as  a  function  of  the  frequency  (u)  for  various  values  of  the  well 
radius  Rc  and  the  dielectric  constant  c.  The  maximum  value  of  the  gain 
(-^  )  in  the  k^-w  diagram  is  shown  in  Fig.  4,  which  determines  the 

maximum  gain  curve  in  the  Rc  -  e  space  in  Fig.  2.  Although  the  chosen 
velocity  spread  (A)  is  1%,  those  for  A:> 3%  are  also  shown  as  a  reference. 

For  A*l%,  the  plot  of  kj  vs.  u  yields  multiple  maxima  for  e  higher  than  a 
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certain  value,  denoted  in  the  figure  with  dotted  lines.  This  multiple 
maxima  phenomenon  is  expected  from  the  considerations  in  Fig.  2, 
corresponding  to  the  case  C.  On  the  other  hand,  for  A=3%,  the  curve 
yield  only  single  peak.  This  can  be  explained  from  the  greater  reduction 
rate  of  the  gain  at  higher  frequencies  for  a  large  spread  [2] ,  [41 ,  [6] . 

That  is,  the  peak  at  higher  frequency  is  quenched  more  rapidly  when  the  spread 
is  large.  The  actual  bandwidth  for  the  same  parameter  in  Fig.  4  is  shown 
in  Fig.  5.  The  bandwidth  is  defined  by  the  full  width  of  the  real  frequency, 
at  which  the  linear  gain  drops  to  exp(-  j)  of  its  maximum  value.  Of  course, 
the  bandwidth  is  normalized  by  its  mean  frequency  u.  This  definition  of  the 
bandwidth  is  somewhat  unconventional,  but  it  serves  the  comparison  purposes. 

For  a  small  spread  (A=l%),  the  bandwidth  curve  in  Fig.  5  yields  interesting 
results.  For  given  wall  radius  CRc) ,  as  the  dielectric  constant  (e)  is 
increased,  the  bandwidth  decreases,  and  then  increases  to  give  a  local  peak, 
followed  by  an  abrupt  decrease.  Although  the  bandwidth  can  be  much  wider  for 
c  lower  than  that  giving  local  maximum,  we  attribute  this  wider  bandwidth 
mainly  to  the  SWM  contribution.  This  is  evident  from  the  bandwidth  curve  for 
A*  3%.  At  e  lower  than  that  giving  local  maxima  for  A=l%,  the  bandwidth  is 
actually  narrower  for  A*3%.  In  view  of  the  more  sensitive  nature  of  the  SWM 
to  the  velocity  spread  compared  to  that  of  the  IWM,  we  therefore  conclude  that 
the  wider  bandwidth  for  low  c  at  A=l%  is  due  to  the  SWM.  Since  our  optimization 
is  for  the  IWM  only  (suppressing  the  SWM),  the  optimized  bandwidth  of  the 
desired  IWM  at  A*H  corresponds  to  the  local  maximum.  The  comparison  of 
Fig.  S  with  Fig.  4  reveals  that  the  maximum  bandwidth  curve  (Fig.  5)  lies 
above  the  maximum  gain  curve  (Fig.  4)  as  expected  from  Fig.  2.  The  optimized 
values  of  the  dielectric  constant  and  the  wall  radius  are  thus  given  by 
(from  Fig.  5)  e*12.6  and  Rc/R®  *  0.77.  These  parameters  are  to  be  compared 
to  c*15.2  and  R£/R®  *  0.63  for  the  wall  clad  configuration. 
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The  dependence  of  the  gain  (and  the  bandwidth)  on  the  bean  center  location 
(R0)  is  illustrated  in  Fig.  6.  Here  the  maximum  gain  is  plotted  against  the 
beam  location  Rq.  The  physical  lower  limit  of  the  beam  location  is  depicted 
by  a  dotted  vertical  line  in  the  figure.  The  maximum  gain  monotonical ly 
decreases  as  the  beam  location  is  increased.  This  dependency  of  the  gain 
is  expected,  since  the  instability  driving  electric  field  tends  to  concentrate 
inside  the  dielectric  center  rod.  Thus,  as  the  beam  is  located  as  close  as 
possible  to  the  edge  of  the  center  rod,  the  gain  is  increased.  At  the  same 
time,  the  SWM  is  also  enhanced  by  locating  the  beam  as  close  as  to  the 
dielectric  material  [2] ,  [4] ,  16] .  This  similar  behavior  of  the  gain  for  the 
SWM  and  IWM  with  respect  to  the  beam  location  is  contrast  to  those  for  the 
wall  clad  case.  There,  since  the  dielectric  is  attached  to  the  outer 
conductor,  the  IWM  gain  is  increased  as  Rfl  increases.  On  the  other  hand, 
the  SWM  gain  is  increased  as  RQ  reduces  (2] ,  [4] ,  [6] .  This  cooperative 

nature  of  the  IWM  and  SWM  in  the  center  rod  geometry  results  in  the  extended 
IWM  region  in  the  frequency  space  as  shown  later.  The  optimized  beam  location 
is  then  given  by  RQ=;  R^  ♦  rL-  The  results  of  Figs.  4-5  are  obtained  at  this 
optimized  beam  location. 

The  process  of  the  bandwidth  optimization  is  further  illustrated  in 
Fig.  7.  Here  the  gain  is  plotted  versus  the  frequency  for  A=l%  and  several 
values  of  the  dielectric  constant.  The  other  parameter  are  optimized  as  shown. 

As  a  reference,  the  case  for  e*l,  that  is  without  the  dielectric,  is  also 
shown  (broken  curve) .  The  maximum  gain  is  achieved  at  the  single  peak  for 
e*12.2,  corresponding  to  a  point  in  the  maximum  gain  curve  in  Fig.  2  (case  B) . 

For  lower  e  (e*11.8)  the  gain  yields  a  single  peak  at  a  lower  gain,  corresponding 
to  case  A  in  Fig.  2.  For  higher  e  (e®12.6,  12.8),  the  beam  mode  intersects 
the  waveguide  mode  at  two  frequencies  (case  C  in  Fig.  2) ,  resulting  in  double 
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maxima  at  those  frequencies.  As  e  increases  from  its  maximum  gain  value 
(ea12.2),  the  gains  at  the  peaks  as  well  as  the  gain  at  the  valley  decrease. 
Therefore,  as  e  increases  from  e=12.2,  the  bandwidth,  now  utilizing  both  peaks, 
increases  until  its  maximum  is  reached  at  e*12.6.  After  the  maximum  bandwidth, 
the  gain  at  the  valley  is  too  low  for  the  bandwidth  to  include  both  peaks 
(e.g.  e=12.8) .  This  explains  why  the  bandwidth  in  Fig.  5  decreases  abruptly 
after  the  maximum  value.  Fig.  7  also  provides  information  on  the  contribution 
of  the  SWM.  Since  the  gain  for  e=l  (without  the  dielectric,  broken  line) 
represents  the  contribution  of  the  SWM  only,  we  can  say  that  any  significant 
difference  from  this  curve  is  due  to  the  IWM.  We  note  from  Fig.  7  that  the 
significant  SWM  contribution  is  for  the  frequency  u>/d>c  >  2.0.  This  is  to 
be  compared  with  u>/u)c  >  1.8  for  the  pure  IWM  and  w/wc  >  1.6  for  the  mixed 
mode  operation  in  the  wall  clad  configuration  (see  Figs.  3  and  4(a)  in 
Ref.  6).  (With  the  wall  clad  configuration,  the  pure  IWM  operation  is 
achieved  with  the  beam  location  close  to  the  wall-clad  dielectric,  and  the 
mixed  mode  operation  with  the  beam  close  to  the  axis,  see  Ref.  6.)  That  is, 
the  SWM  begins  to  contribute  significantly  to  the  bandwidth  at  considerably 
higher  frequency  with  the  center  rod  configuration  than  with  the  wall  clad 
configuration,  thereby  enhancing  the  frequency  range  of  the  IWM  for  the 
center  rod  configuration.  In  view  of  the  difficulties  associated  with  the 
SWM,  therefore,  the  center  rod  configuration  is  superior  to  the  wall  clad 
geometry. 

The  foregoing  optimization  process  with  respect  to  the  bandwidth  for  1% 
of  the  axial  velocity  spread  and  the  beam  parameters  given  in  Eq.  (9)  can  be 
summarized  as  following: 

R ,/R_  *  0.20,  e  -  12.6,  R  /R°  *  0.77 
d  c  c  c 

R0  «  Rd  +  rL  -  0.35  Rc.  (11) 
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The  optimized  bandwidth  at  A*l%  is  68%  (Fig.  7),  of  which  51%  is  due  to  the 

INM,  and  17%  to  the  SIM.  This  proportion  is  obtained  from  Fig.  7, 

attributing  the  portion  for  <n/uc  >_  2.0  to  the  SNM.  This  bandwidth  can 

be  compared  to  46%  (all  INM)  for  the  pure  IWM  and  90%  (12%  INM  ♦  78%  SUM) 

for  the  mixed  mode  with  the  wall  clad  configuration  [6].  Moreover,  the 

mean  frequency  u  for  the  center  rod  gyroton  is  u/u  ■  1.72,  higher  than  1.53 

c 

for  the  pure  INM  and  less  than  2.04  for  the  mixed  mode  with  the  wall  clad 
geometry  [6].  If  we  assume  that  the  entire  SWM  contribution  is  non-usable, 
then  the  bandwidth  of  tK$  center  rod  gyrotron  is  slightly  wider  than  that  of 
the  pure  INM  and  much  wj^er  than  that  of  the  mixed  mode  of  the  wall  clad 
gyrotron. 
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PERTURBED  FIELDS 


In  previous  section  we  have  seen  that  the  frequency  range  of  the  IWM 

for  the  center  rod  configuration  is  wider  than  that  for  the  wall  clad 

geometry.  Since  the  IWM  results  from  the  beam-waveguide  mode  coupling,  its 

perturbed  field  profile  is  very  similar  to  that  of  the  beam- free  waveguide, 

that  is,  of  Bessel  function  type.  On  the  other  hand,  the  field  profile  due 

to  the  SWM  is  highly  localized  near  the  beam  location  [2] ,  [4] ,  [6] . 

Therefore,  we  can  determine  the  nature  of  the  instability  by  plotting  the  field 

profiles.  The  instability  driving  fields  E0J  and  BrJ  are  given  in  Eqs.  (3) 

and  (4).  Since  the  instability  occurs  near  the  beam  mode,  the  field  profiles 

are  computed  at  k*k„  =  (w-u  /y)/c8,  for  given  u>.  The  quantity  k_  is  the 

wavenumber  corresponding  to  u-Ug  (Eq.  (8)).  In  the  IWM  frequency  range,  it  is 
2  2  2 

shown  that  (u  /c  -  kR)  <  0,  hence  the  Bessel  functions  J  and  N  with  arguments 

x's  (Eq.  (7))  now  become  the  modified  Bessel  functions  I  and  K. 

In  Fig.  8,  the  field  profiles  of  E01  and  Brl  at  a  particular  frequency 

us1.7  u  for  several  values  of  the  dielectric  constant  are  shown.  That  is, 
c 

the  profiles  are  drawn  from  Fig.  7  at  that  frequency.  The  fields  are  normalized 
such  that  J"  c  Ej  ♦  Bj)  rdr  *  R^.  Here  *  e(1)  for  r  <  Rd  (r  >  Rd) .  The 
chosen  frequency  <u*1.7  corresponds  approximately  to  the  transition  point 
between  the  IWM  and  SWM  (see  Fig.  7) .  We  observe  in  Fig.  8  that  the  profile 
for  e*l  (without  the  dielectric)  is  highly  localized  at  the  beam  location  (Rq) 
with  a  negligible  field  amplitude  in  the  space  (r  <  R^)  where  the  dielectric 
rod  is  supposed  to  be.  This  is,  of  course,  expected  since  the  instability 
for  e*l  is  solely  due  to  the  SWM.  Without  the  dielectric  (c*l),  the  waveguide 
mode  is  a  fast  wave  (o>G  >  ck) ,  and  the  slow  wave  IWM  is  absent.  As  c  increases 
the  SWM  contribution  is  decreased,  as  indicated  by  the  decreasing  peaks  of  the 
field  profile  at  the  beam  location.  The  reduced  SWM  contribution  with  the 
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increased  e,  in  turn,  indicates  the  more  enhanced  IWM  contribution  (see  Fig.  2). 
For  example,  at  e*12.8,  the  field  profile  is  almost  identical  to  that  of  the 
waveguide,  indicating  nearly  pure  IWM  contribution  to  the  instability.  At 
u/o>c  *  1.7,  the  significant  SWM  contribution  is  for  e  <  11.8,  as  can  be 
confined  from  Fig.  7.  In  terms  of  difficulties  associated  with  the 
microwave  excitation  and  collection,  the  localization  of  the  field  strength 
due  to  the  SMM  at  the  beam  location  can  be  a  nuisance.  However,  these 
difficulties  can  be  eliminated  by  optimizing  the  parameter  e,  i.e.  e=12.6 
in  Fig.  8. 

The  frequency  dependenc  of  the  field  profile  is  shown  in  Fig.  9  for 
the  optimized  parameters  given  in  Eq.  (11).  As  the  frequency  increases,  the 
peak  at  the  beam  location  (RQ)  is  more  pronounced.  This  indicates  that  the 
contribution  of  the  SWM  is  nearly  negligible  for  u/u c  <  l-8  and  is  significant 
for  u/uc  >  2.0,  confirming  the  similar  conclusion  from  Fig.  7.  This  profile 
can  be  compared  with  that  for  the  mixed  mode  with  the  wall  clad  geometry 
(Fig.  6  in  Ref.  6),  where  the  SWM  is  significant  for  w/o>c  >  1.6.  We  thus 
again  find  that  the  center  rod  configuration  extends  the  IWM  region  further 


than  the  wall  clad  gyrotron. 
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CONCLUSION 

We  have  investigated  the  wide  band  capability  of  the  gyrotron  with  a 
dielectric  material  used  as  a  center  rod.  After  deriving  the  dispersion 
relation  for  the  azimuthally  symnetric,  TE  perturbations,  we  have  found  the 
optimization  conditions  on  the  physical  parameters  for  a  wide  bandwidth  at  a 
small  axial  velocity  spread  (A=l%).  ' 

The  results  of  the  optimization  processes  and  the  comparison  with  the 
wall  clad  configuration  can  be  summarized  in  Table  1.  All  values  are 
obtained  with  the  beam  parameters  in  Eq.  (9).  The  center  rod  configuration 
allows  larger  wall  dimension,  lower  dielectric  constant,  and  thicker  center 
rod  compared  to  the  wall  clad  configuration.  The  gain  for  the  center  rod 
configuration  is  slightly  lower  than  that  of  the  pure  IWM  for  the  wall  clad 
configuration,  but  higher  than  that  of  the  mixed  mode  for  the  wall  clad 
configuration.  The  bandwidth  is  wider  than  the  pure  IWM,  but  narrower  than 
the  mixed  mode.  However,  in  terms  of  difficulties  associated  with  the  SWM, 
the  center  rod  configuration  proves  to  be  superior  in  the  bandwidth  to  both 
the  pure  IWM  and  the  mixed  mode.  This  advantage  comes  from  the  extended  IWM 
frequency  range  in  the  center  rod  configuration. 
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Figure 


1  Cross  section  of  a  gyrotron  with  a  dielectric  center  rod 
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NOTE:  The  dispersion  curves  for  the  parameters  represented  by 

A,  B  and  C  are  shown  in  the  upper  corner^  The  maximum  bandwidth 
curve  (broken)  is  located  above  the  maximum  gain  curve  in  Rc~c 
space. 


Figure  2  Schematic  diagram  of  the  mode  characteristics  in  the  space  of 
the  conducting  wall  radius  (Rc)  and  the  dielectric  constant  (c) 
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NOTE:  The  values  of  Rc  and  e  that  make  the  beam  mode  graze  the  waveguide 
mode  (upper  corner)  are  shown  for  several  values  of  the  thlckitaod 
parameter  (R^/Rc> . 


Figure  3  Grazing  conditions  on  the  conducting  wall  radius  (Rc,-)  and  the 


dielectric  constant  (c,— -) 
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NOTE:  For  A«l%,  the  gain  exhibits  double  aaxiaa 
at  high  e  denoted  by  broken  lines. 


Figure  4 


Plots  of  the  aaxinum  gain  (-k^  )  vs.  the  dielectric  constant  (c) 
for  several  values  of  the  conducting  wall  radius  (Rc),  and  at  two 
different  axial  velocity  spreads  (A) 


29 


1 


— 100  k,  c/co 


NCWC  TM  «1-at7 


NOTE:  The  broken  vertical  line  corresponds  to  the  physical  lower 
limit  of  Rq,  where  the  inner  edge  of  the  beam  touches  the 
dielectric  rod. 


Figure  6  Dependence  of  the  maxima  gain  (-k^**)  on  the  beam  center 
location  (Rq) 
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Figure  7 


Plots  of  the  gain  (-k^  vs.  the  frequency  (u)  for  several  values 
of  the  dielectric  constant  (c).  Other  optiaized  paraaeters  are 
as  shown.  The  broken  line  (e  ■  1.0)  corresponds  to  the  absence 


of  the  dielectric  rod 
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Figure  i  Perturbed  field  profiles  for  several  values  of  the  dielectric 

constant  (c)  for  u/«c  "1.7  and  parameters  otherwise  identical  to 
those  in  Fig.  7 
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Table  1  Comparison  of  the  bandwidth-optimized 
dielectric  gyrotron 


Dielectric 

Used  As 

Center  Rod 

Wall  Clad 

Pure  IWM 

Mixed  Mode 

£ 

12.6 

15.2 

VRc 

0.77 

0.63 

Thickness  Ratio 

R./R  =0.20 

a  C 

VRc  =  0>85 

Beam  Location 

R0  i  Rd  *  rL 

RQ/Rc  =0.38 

-100  k“x  c/uc 

4.6 

5.6 

3.8 

Significant  SWM 

Contribution 

o 

<M 

A  1 

O 

3 

3 

<*»/  u>c  >  1.8 

<*>/  (i)  >1.6 

c 

Mean  Frequency 

w/u  ■  1.72 
c 

w/w„  *  1.53 
c 

fc/w_  »  2.04 
c 

Bandwidth 

(INM  +  SWM) 

68% 

(51  ♦  17) 

46% 

(46  ♦  0) 

90% 

(12  ♦  78) 

-100  k®45'  c/w 

2.7 

3.9 

1.8 
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